arXiv:hep-th/9812107vl 14 Dec 1998 


Quasi-relativistic center-of-mass variables 
in applications 


R.Gaida , V.Tretyak, Yu.Yaremko 


Institute for Condensed Matter Physics of the 
National Academy of Sciences of Ukraine 
1 Svientsitskii St., UA-290011 Lviv, Ukraine 


Collective center-of-mass variables are introduced in the Lagrangian formal¬ 
ism of the relativistic classical mechanics of directly interacting particles. It is 
shown that the transition to the Hamiltonian formalism leads to the Bakamjian- 
Thomas model. The quantum-mechanical system consisting of two spinless par¬ 
ticles is investigated. Quasi-relativistic corrections to the discrete energy spec¬ 
trum are calculated for some Coulomb-like interactions having field theoretical 
analogues. 
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1. Introduction 

There are two types of classical relativistic direct interaction theories (RDIT), 
namely, manifestly Poincare-invariant four-dimensional formulations and three-di¬ 
mensional approaches [ |l| . Applying each of them to the description of real physi¬ 
cal Af-particle systems we come across the problem of space-time interpretation of 
some formal constructions. In four-dimensional approaches difficulties are caused 
by the presence of redundant variables specific to many-time formalisms and in 
three-dimensional ones they are caused by the so-called no-interaction theorem. In 
the relativistic Hamiltonian mechanics, this theorem forbids identification (in the 
interaction zone) of canonical variables q a l with physical (covariant) position coor¬ 
dinates Xa (a — 1 , N‘ i — 1,2,3) which transform as a space part of the events 
in the Minkowski space-time [ |2[. In the Lagrangian formulation [ a a. u a the 
no-interaction theorem asserts that the interaction Lagrangian depends inevitably 
on all higher time derivatives x a l = d s x a l /dt s , s = 1 , 2 ..., of the physical parti¬ 
cle coordinates. This theorem is essentially related to the use of individual particle 
variables (PV) in both the Hamiltonian and Lagrangian approaches. 

In papers [BSS11 the following approach to the construction of RDIT was 
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proposed. The starting point is the 3-dimensional Poincare-invariant Lagrangian 
formalism ib a with using the physical coordinates x a l and with allowing the 
connection to field theories via the Fokker-type formalism [ ||(|]. These aspects of 
the Lagrangian description along with the possibility to use the Not her theorem 
for obtaining the conservation laws [ (4[ are the main advantages of this approach 
from the point of view of the physical interpretation of the theory. The next step is 
a transition to the canonical Hamiltonian formalism admitting the quantization of 
the theory [111- 

Another aspect of RDIT is the problem of how to separate the motion of the 
system as a whole from the internal motion of its constituents. This problem can 
be solved in a manner similar to that of the non-relativistic mechanics case by 
introducing center-of-mass variables (CMV). Relativistic generalization of the non- 
relativistic notion of the center-of-mass is ambiguous. The most frequently used 
are, firstly, the center-of-inertia (covariant center-of-mass) whose coordinates R l are 
transformed by the Poincare group as 3-coordinates of the Minkowski-space points, 
and, secondly, the center-of-spin (or canonical center-of-mass) with the canonical 
coordinates Q l [||- External variables Q l and P % describing the motion of a system 
as a whole, along with the corresponding set ( qb\Pbi ), b — 1,N — 1 of 6(iV — 1) 
internal variables are often used in the relativistic Hamiltonian mechanics [ [T0|-[ 

PI. 


In the relativistic Lagrangian formalism it is natural to choose the coordinates of 
the center-of-inertia R as external position variables and to formulate the problem of 
constructing the relativistic Lagrangian mechanics of an iV-particle system in terms 
of CMV. This set also contains 3(N — 1) internal variables /V [0-[0- The internal 
variables /V may be introduced by postulating their transformation properties with 
respect to the Poincare group. The corresponding relations which are written down in 
section 2 constitute the prerequisites which enable us to find a connection between 
CMV and PV and, accordingly, between the relativistic Lagrangians in terms of 
these two classes of variables. An important feature of CMV is the existence of the 
Poincare-invariant interaction Lagrangians which depend only on a finite number of 
time derivatives (section 3). The transition to the Hamiltonian formalism (section 
4) leads to the well-known Bakamjian-Thomas model supplemented by the relations 
between the canonical variables and the physical particle positions. In sections 5 and 
6 the first-order quasi-relativistic approximation (up to c -2 terms) is considered from 
both the classical and the quantum mechanical viewpoints. Relativistic corrections 
to the energy spectrum of a two-particle system with Coulomb-like interactions are 
calculated as an illustration. 


2. Realization of the Poincare group V(l,3) in terms of the La¬ 
grangian CMV 

Since the external variables R l describe the motion of an iV-particle system as a 
whole, we require their transformation properties with respect to the Poincare group 
to be the same as those of the physical coordinates x a l of a single point particle. In 
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the instant form of dynamics [ |j, [1] they are [ |L|, [17 


TiRi = R{ , 

K-iRj 


ViRj &ij i JiRj £ijkRk • 

—tSij + c RiRj. 


( 2 . 1 ) 

( 2 . 2 ) 


Here H := D — d/dt, Vi , J7), and /Q are generators of time and space translations, 
space rotations, and boosts, respectively, and D is the total time derivation operator. 
As usual, the dot denotes the time derivation (e.g., Rj = DRi are the components 
of the external velocity R). We postulate further that similar to the case of non- 
relativistic mechanics the internal CMV p b are the components of the translationally 
invariant 3-vectors p b (b — 1, N — 1) which do not depend explicitly on time: 


R-Pbi — Pbi i ViPbj 0 , \JiPbj — £ijkPbk ■ (2.3) 

The transformation of these variables with respect to boosts is determined by the 
Poincare group structure. It can be shown that the commutation relations of the 
Poincare algebra ^4.7^(1, 3), along with equations ( p.l| )~( p3[) , lead to the following 
conditions [ 

R'iPbj C ^RiPbj T RiPbj Aj (Rpj f ). (2.4) 


14, 


The expressions for generators X a = {7i, Vi, /Q} of the Poincare group in terms 
of CMV may be obtained from equations (|2.1|) - (|2.4|) immediately [ [T| . 


3. Poincare-invariance of the Lagrangian description 

The evolution of a particle system is described by the solutions of the Euler- 
Lagrange-Ostrogradsky equations 

°° d I 

E<- b )‘-^ = o. KmW.rtT (3.1) 

s=0 o Z„ 

following from the variational principle 5 J Ldt = 0. The conditions of Poincare 
invariance of equations (|3.1|) may be expressed by the system of first-order differential 
equations for the Lagrange function: 

X a L = DQ a , (3.2) 


where the auxiliary functions have to satisfy the compatibility conditions: 

XaVlf) — , (3.3) 

with c Q p 7 being the structure constant tensor of the Poincare group V(l, 3). For this 
group the relations (|3.3|) admit the following set of functions 


H-h — L , Q-p t — 0 , — 0 , QjCi — c 2 RiL 


(3.4) 
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In which case the conditions (|3.2|) of Poincare invariance take on the form: 


V t L = 0, 


JiL = 0, (H-D)L = -— = 0, 
(/C i — c~ 2 RiD)L — c~ 2 RiL = 0 . 


(3.5) 


The structure of equations is quite different from the structure of the 

Poincare invariance conditions in terms of PV x a l | []]. [4|j. Whereas the interaction 
Lagrangian as a function of the particle variables is defined on J°° 7 r, equations ( |3.5| ) 
admit Lagrangians that depend on CMV and their time derivatives of finite orders. 
The simplest structure of L is given by 

L{ R, R, p b , p b ) = T~ 1 F(a b , (3 b , 7 ), b = 1 , N - 1 ; (3.6) 

here T ” 1 = \/l — c _ 2 R 2 and F is an arbitrary function of scalar products of the 
vectors 


ct b = p b T c ^[Rp 


/3 b = Fct b — c 1 r 2 [RcK 


7 = DR + c _ 1 r 2 [RR] + c _ 2 r(RR)R. 


(3.7) 

(3.8) 


Contrary to ct b and f3 b , the 7 dependence of L is not necessary. The Lagrangian ( |3.6| ) 
depending on a finite number of CMV can be made equivalent to the Lagrangian 
defined on an infinite-dimensional set of PV x a 1 , only if the CMV are functions of 
all these PV. The set of equations (|2T| ) -(|2~^ ) with generators H, Vi, Ji, and /Q, 
determined in terms of PV (IB- may be regarded as the basis for determining 
these functions [ O . 

Using the Nother theorem for Lagrangians with higher derivatives [|-[| 


we 


found ten integrals of motion (E, P, J, K) corresponding to the Poincare-invariance 
of the theory, based on Lagrangian ( |3.6| ). From their analysis we obtain [ [L4], [L7]] the 
following equation of external motion 


MR + c _ 2 [R, s] = K 0 


(3.9) 


in the rest frame of reference (E = E 0 = Me 2 , P = 0, J = s, K = K 0 )). It should 
be noted that the equation 


MR + c -2 [R, s] = 0, 


(3.10) 


being the consequence of ( |3.9| ), was obtained in [ [L 8 |] for classical particles with a 
spin; it is a three-dimensional form of the well-known Mathisson equation. Spin s is 
equal to the internal angular momentum of the system. 

Solving ( 11.91) in a coordinate system with the origin at point R 0 = K 0 /M, where 
K 0 is the value of integral K in the rest frame and the CU axis is directed along 
spin s, we see that point R(t) describes a circular motion in a plane orthogonal to 
the spin: 


R x (t) = Asin(u;ot — a), R y (t ) = Ylcos(u;ot — a), R z = 0 . 


(3.11) 
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Here, frequency uq = Mc 2 /s; A and a are constants of the integration. Taking 
into account the usual non-relativistic limit procedure, in accordance with which 
the center of mass is at rest as measured in the rest frame, we must set A — 0. 
Indeed, as c~ 2 —> 0 (cno —»■ oo), the limits of the trigonometric functions in equation 
( [3.11D do not exist. This means that in the set of solutions of equations (pT|) for 
the Lagrangian ( |3.6|) we keep only those which satisfy the standard free-particle 
equation R = 0. 


4. The Hamiltonian description 


There are several ways to obtain the Hamiltonian formulation corresponding to 
the Lagrangians with higher derivatives and with supplementary conditions, like 
R = o [ g [11)1 PH] • Here we use an approach close to that of [ [l9| . Let us perform a 
change of the variables (R, p b ) i—► (Q, q&), where 


Q = R 


q b = r p b 


[P-s] 


M(E + Me 2 ) ’ 

c ~ 2 r 2 ■ . 

-R(R P b ) 


r + i 


(4.1) 

(4.2) 


Relation (|4.1| ) is a definition of the canonical center-of-mass [ |j. In terms of new 
variables, the Lagrangian ( ]'i.G| ) has the form: 

L(R, R, p b , p b ) = L(Q, 0, q b , q b ) + U = L{ Q, q b , q b ) + U , (4.3) 

where L' contains only the terms which do not contribute to the equations of mo¬ 
tion if the external equation of motion R = 0 is satisfied. Therefore, the term U 
may be omitted in equation (|4.3|) and, thus, the transformation to the Hamiltonian 
description may be performed in the usual manner on the basis of Lagrangian L. 
Using the Legendre transformation we construct Hamiltonian H corresponding to 
Lagrangian L : 

H = a/c 2 P 2 + M 2 (q b , 7T b )c 4 , (4.4) 

where 


dL 

9Q 


7U = 


dL 

dq b 


N —1 o j-i 

Me 2 = £ _ F 


6=1 


•S0 b 


(4.5) 


R=0 


We note that the external momentum variable R is a motion integral (the total 
momentum of the system). 

Expression (|4.4| ) for H reproduces the particle-system Hamiltonian of the Ba- 
kamjian-Thomas model m . The advantage of our way of constructing this model 
is the possibility to express the internal Hamiltonian M(q b , n b ) in terms of the 
interaction Lagrangian which, in turn, may be related to some field models [ |S). On 
the other hand, in our approach the canonical variables Q.P.qj,, and ir b have a 
definite relation to the Lagrangian CMV R, p b and their time derivatives R, p b and, 
therefore, to the individual particle variables x a , at least in the approximations in 
the parameter c~ 2 . 
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5. Quasi-relativistic approximation 


Let us consider the main steps of our approach in the first-order approximation in 
c -2 (i.e. quasi-relativistic or post-Newtonian approximation). We use the following 
sufficiently general form of a quasi-relativistic two-particle Lagrangian with static 
non-relativistic potential U (r) I00i : 

a= 1 ' ' 

+ ((rv!)(rv 2 ) + 25(rv) 2 ) + 2«i(r)^ , (5.1) 


where r = xi — x 2 , v = Vi — v 2 , A and B are arbitrary constants and U\{r) is an 
arbitrary function. For interactions which correspond to the linear field theories, we 
have u\{r) = 0 and constants A and B depend on the tensor structure of interaction [ 
|1], [2TJ] . For example, we have A = B = 0 for a vector interaction, and A = — 1 , B = 0 
for a scalar interaction. 

Further we consider a specific case of interactions which correspond to the non- 
relativistic long-range potential 

U(r) = ~, (5.2) 

r 

where a is a constant of the interaction. We suppose that the Galileo-invariant 
function U\{r) has the form: 


/ \ ^ 1 
“iM = 72 • 


(5.3) 


where aq is the constant. Let us note five significant cases covered by the interaction 
described in a quasi-relativistic approximation by the Lagrangian of type ( |5.1|) where 
functions (|5.2|) and ( |5.3|) enter: 


• a = eie 2 , A = B = 0, and aq = 0: the Darwin Lagrangian which describes an 
electromagnetic interaction of two point charges; 

• a = -Gm 1 m 2 , 2 A — 3, B — 0, and 2aq = G 2 mim 2 m(m = m\ + m 2 ) : the 
Einstein-Infeld-Hoffmann Lagrangian which describes a gravitational interac¬ 
tion of two point masses; 


• a = eie 2 — Gmi7n,2,2A = — 3Gm 1 m 2 /cx, B = 0, and 2aq = G(Gmim 2 m + 
e 2 m 2 + e\mi — 2e\ e 2 m): a gravitational interaction of charged bodies (the 
Bazanski Lagrangian [0)i 

• a = -Gm 1 m 2 , 2 A = 1 + 2y, B = 0, and 2aq = ( 2(5 — 1 )G 2 m 1 ?n 2 m: a gravita¬ 
tional interaction in the PPN-formalism [ ^3|] which classifies suitable theories 
of gravitation (i.e. those which do not contradict the four classical experi¬ 
ments, satisfy the conditions of the Lorentz-invariance, and admit IV-particle 
Lagrangian functions) in the first-order approximation in c” 2 ; 
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• a = —Gmim 2 ,A = —2 B = 1, and ay = 0: the Lagrangian of the quasi- 
relativistic approximation of the Whitehead-Shield theory of gravitation [ 
where the massless tensor held of the second rank as an intermediate held is 
assumed. 


The correspondence between CMV and PV for a two-particle system in the first- 


order qnasi-relativistic approximation has the following form [ [17 


x„ = R + (-l)*^p+i((-l)^(iR 2 p-(Rp)R]- 
m c 2 m \ 2 


m 


^(Rp)p + 


m- 


m i — m 2 
2m? 


(/'P J - U(pj) p \, 


(5.4) 


where a,b = 1,2, 6 ^ a,m = m\ + m 2 , p = mim 2 /m. Using these relations in the 
Lagrangian ( |5.1|) leads to the following two-particle quasi-relativistic Lagrangian in 
terms of CMV: 


L = 


mil 2 upi TT , . mR 1 up 

— + ~TT- u (p) + + 


8c 2 


P /Y> a\2 , P 


2 c 2 


(Rp ) 2 + ^ R[p[Rp]] + 


8c 2 

1 

2d 2 


1 - - V T^ 2 P 2 - 

m / 4 cr 


R 2 - 


P_ 

m 


+ 2A) p 2 ) U(p)+ 




(5.5) 


This expression agrees with the corresponding approximation of the Lagrangian 

ra¬ 


in the quasi-relativistic approximation the canonical external variable (|hl|) be¬ 


comes 


Q 


R + 
R + 


[P.sl 


M(E + Me 2 ) 


P 


2 me 2 - 


R[pp] 


R=0 


(5.6) 


The total momentum 

P = P|r.o = fa (5.7) 

is canonically conjugated to Q. Quasi-relativistic inner canonical variables are de¬ 
fined by the relations 


q 


7T 



pc 2 dp ’ 

_ ^dA 
c 2 dp 


(5.8) 

(5.9) 
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where \(p, pi) is an arbitrary Galileo-invariant function. Choosing this function in 
the form 


Hp.p) = (pp) I (i - - 

m 


+ i( 1 -s- 4fl ) ,7W )- 


(5.10) 


we obtain the Hamiltonian function depending on the square of the inner momentum 
7r and therefore being convenient for quantization: 


H = 


P 2 ( 7T 2 

V + U(q) 


2m 2 m 2 c 2 \2/i 


P 4 7T 2 

+ 


8 m 3 c 2 2 p(q) 


v(q) . 


Here 


p{q) ijl ' 4/i 2 c 2 1 2 ( 1+ m + 4A 4B ) ( X m ) Q dq 


and 


%) = U( q ) + T \ uM - ^ (l - T _ 4fi) C(5)9^ 


dU 

dq 


(5.11) 


(5.12) 


(5.13) 


6. Quasi-relativistic Schrodinger equation 


Starting with the Hamiltonian function (|5.11|) we write down the stationary 
Schrodinger equation: 

(H — E)(Q, q | E.P) = 0 . (6.1) 

The total wave function (Q,q|_E',P) of the system is equal to the product T(Q)-?/>(q) 
of the external wave function T(Q) and the inner one tZ’(q). We consider the states 
with sharply defined total momentum P. so the external wave function is then 
a plane wave, d'(Q) = exp(ikQ), where the wave vector k = P/h. Due to the 
separation of the external and internal variables in equation ( p.l| ), we obtain the 
following wave equation: 

(h-Sm q) = 0, (6.2) 


where 


_ P 2 EP 2 
2 m 2m? c 2 


8m 3 c 2 ' ^ 6 ' 3 ^ 

Here E and P are eigenvalues of the non-relativistic Hamiltonian and total momen¬ 
tum operators, respectively. 

The internal Hamiltonian operator is defined by the formula 


ft = -A' + v(q) , 


(6.4) 


where the first term is some Hermitian operator which corresponds to the classical 
function vr 2 / / n(g) = a(q)7T 2 . We use the following sufficiently general quantization 
rule: 

X = a(q) ft 2 + 5 lj [jfi, a(q)]7tj + Ad 2 (a(g)) , (6.5) 
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where the inner momentum operator is chosen in the usual form: 

d 


iti = —ih 


dq i 


r = -n 2 a . 


( 6 . 6 ) 


This rule preserves the rotational invariance of the inner Hamiltonian and depends 
on one parameter A G M. The value A = 1/2 corresponds to the symmetrization 
( 7 f 2 a(q) + a(q) 7t 2 )/2. 

To illustrate our approach we consider the Coulomb-like potential (see equations 
( |5.2D and (|5.3|) ). In this case the functions a(q) := fJL ~ l (q ) and v(q) take the form: 


a(g) = 


v(q) = 


/i 


a 


Q 


1 + -C,- 
lic z q 

1 + \c 2 - 

/icr q 


Ci = Vik +2A - 2B ’ {ej) 

= + (6.8) 

4 4m cr 


Here /i = m\m 2 /m and the parameters A, B, cti, and a are defined by the type of 
interaction (see section 5). 

After some calculations we arrive at the following Hamiltonian operator: 


h = 


h 2 


1 . a _ q d ATTaX^ . . 

A -|——Ci— 7 ;—I-Tr^rC7i(5(q) ) + v(q ), 


Kq) 


fi 2 c 2 q 3 dq fi 2 c 2 


(6.9) 


where 7r denotes a well-known transcendental number instead of the value of the 
inner momentum. 

We shall study the discrete spectrum of H, when a < 0. The use of methods of the 
perturbation theory (e.g., see [|25|]) allows us to obtain quasi-relativistic corrections 
to the non-relativistic Bohr expression 


£i 0) = 


/iCX 2 1 

2 h 2 n 2 ’ 


( 6 . 10 ) 


where n is the main quantum number. Let us rewrite the Hamiltonian operator ( |6.9| ) 
in the form 


h = 




Me) 




C-> — C\ a 2 fi 2 a „ ( q d 


/j,c 2 q 2 


+ 


where 


2 ^ Cl (ysi +4,rA ' S(q) 


h 2 


U °) = ——A + - 
2/i q 


( 6 . 11 ) 


( 6 . 12 ) 


is a non-relativistic energy operator possessing spectrum ( |6.10| ). According to the 
stationary perturbation theory [ |25|, the quasi-relativistic term is equal to the 
matrix element 

(6.13) 


£nl = I = (V)o, 


,,*( 0 ) 


( 0 ) 
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where q) are the eigenfunctions of operator ( |6.12| ) and 


v = T Cl “ft<»> + Az^I + ^_ Cl 


/iC z 


q 


/iCY q- 


2 jj?c 


+ 47rA5(q) 

q A <Aq 


(6.14) 


is a quasi-relativistic perturbation of the non-relativistic energy operator ( |6.12D . 
While calculating (|6.13|) we take into consideration that //’^(q) are the eigenfunc¬ 
tions of hS 0) . We also use the average values of the following functions of r = q [ 

E3: 


<r- 2 )o = 

n z n z 

(5(r)) 0 = IV'SL 


/i 2 ct 2 


h 4 n 3 (2/ + 1) ’ 


2 ^l«| 3 


irh 6 n 3 

One can easily check that the average value is 


oi ■ 


d 3 #X } (q)^^^SL(q) = 2tt|^1(0)| 2 . 


(6.15) 


(6.16) 


Therefore, the quasi-relativistic correction to the energy spectrum is given by 

(6.17) 


c(i) ^ 1 (Cl 4 (C!-C 2 ) , ox lN 2C lr 

£ ni = oTToii ft -wo-( 2A - !)— A 


2 h 4 c 2 n 3 \ n 


21 + 1 


n° 


0 / 


We note a dependence of the obtained spectrum on the quantization parameter A. 

For the case of an electromagnetic interaction, when the parameters A = B = 0, 
a>i = 0 (see section 5), spectrum ( |6.17| ) agrees with the results of the conventional 
quantum electrodynamics provided that A = 1/2. Let us consider a hydrogen-like 
atom. Then the constant of interaction a = —Ze 2 , where Ze is the charge of the 
nucleus. Since the mass of the nucleus is much greater than the rest mass of electron 
m e , we have /j —» m e and, therefore, the ratio fx/m —> 0. Here m is the total 
mass of the atom. Introducing these parameters into equations ( |6.7| ) and (|6.8|) we 


obtain C\ = 3/4 and C 2 = 1/4. According to equation (|6.17|) , the quasi-relativistic 
correction is equal to 


c(i) _ 

L 'nl ~ 


m e {Ze 2 ) 4 1 


h 4 c 2 


rr 


n 


4 V 8 21 + 1 


(6.18) 


It is in good agreement with the quasi-relativistic correction to the hydrogen-like 
spectrum which follows from the Dirac equation (see [ |25], p.216]). 

Now let us consider parapositronium. The mass of an electron is equal to the 
mass of a positron and, therefore, fi = m e /2 and the ratio /i/m = 1/4. Whence the 
constants C\ = 11/16 and C 2 = 3/16. Putting these parameters into equations ( |6.17| ) 
we re-obtain an expression for the quasi-relativistic splitting of the parapositronium 
spectrum which is conventionally found by using the Breit equation [ |26|: 


cd) _ 

°nl — 


771 , 


: (e 2 ) 4 1 ( 11 


2 h 4 c 2 n 3 V 32n 21 + 1 


(6.19) 
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7. Conclusions 

The transformation of particle variables into center-of-mass ones for a class of 
two-particle quasi-relativistic Lagrangians with interactions having field-theoretical 
analogues has been carried out. As a result, the motion of a system as a whole was 
separated from its inner motion. A set of canonical variables in terms of which the 
Hamiltonian has a quadratic dependence on the inner momentum variable is found. 
The problem is quantized. The quasi-relativistic splitting of the discrete energy 
spectrum is obtained. As expected, the degeneracy in the orbital quantum number 
is eliminated. 

In the specific case of an electromagnetic interaction between particles, the quasi- 
relativistic corrections to the non-relativistic discrete spectrum are in agreement with 
the corresponding approximation of the spectrum of a hydrogen-like atom obtained 
by using the Dirac equation (modulus of the spin of the electron which was not 
taken into account). For parapositronium (where the summary spin is equal to zero) 
the spectrum obtained within the frame of the CMV-formalism coincides with the 
one derived by using the Breit equation. 
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